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INTRODUCTION 


Many  very  large  16-inch  inner  diameter  ( ID)-jacketed  thick-walled  pressure 
vessels  manufactured  about  1945  have  recently  been  returned  to  service.  These 
vessels  are  of  a  jacketed  design  as  shown  schematically  in  Figure  1.  Long,  but 
shallow,  longitudinal  10  cracks  have  appeared  in  some  of  these  vessels  during 
repeated  pressurization.  The  purpose  of  this  report  is  to  establish  the 
tolerance  of  the  vessels  to  the  presence  of  these  cracks. 

PROCEDURE 

To  establish  the  damage  tolerance  of  any  structure,  first  the  material's 
fracture  properties  must  be  determined.  Second,  the  stress  intensity  factors 
that  are  applied  to  the  structure  during  load  application  must  be  calculated. 

And  finally,  using  both  of  these  results  the  critical  flaw  size  and  the  number 
of  loadings  that  will  grow  the  crack  to  that  size  can  be  predicted. 

The  fracture  mechanics  properties  required  are  the  fracture  toughness  of 
the  material  (Kjc)  and  the  fatigue  crack  growth  rate  as  a  function  of  the 
applied  stress  intensity  factor  (dl/dN  versus  AK).  The  fracture  toughness  prop¬ 
erties  were  measured  from  the  largest  samples  that  could  be  machined  from  the 
liner  material  (the  liner  is  the  innermost  segment  of  the  vessel  and  is  the 
structural  component  that  cracked).  The  specimen  design  was  the  compact 
tension-type  sample  with  a  through-thickness  dimension  of  1.5  inches 
[C(T)-1.5T].  The  crack  propagation  data  were  alsa  obtained  from  the  liner 
material.  The  specimen  design  for  this  was  a  pure  bending  sample.  This  speci¬ 
men  could  be  loaded  so  that  negative  loading  ratios  could  be  applied  which  would 
simulate  the  effects  of  compressive  residual  stress  due  to  the  interference  fit 
between  the  liner  and  the  outer  jackets  of  the  vessel.  Crack  growth  rate  data 


were  obtained  at  two  loading  ratios,  one  where  the  loads  went  from  zero  to  ten¬ 
sion  loading  and  the  other  for  fully-reversed  loading. 

As  previously  stated,  in  order  to  predict  critical  flaw  sizes  and  determine 
remaining  life,  the  stress  intensity  factor  solution  for  the  structure  under 
load  must  be  determined.  For  simple  cases,  this  can  either  be  estimated  or 
found  in  the  literature.  The  K  solution  for  this  specific  geometry  was  devel¬ 
oped  using  finite  element  analysis,  since  there  were  no  solutions  in  the  litera¬ 
ture  for  the  specific  geometry  of  the  16-inch  vessel. 

The  finite  element  routine  used  was  the  ABAQUS  (ref  1)  general  purpose 
program.  This  code  was  developed  to  analyze  nonlinear  problems,  but  it  also  has 
an  excellent  linear  analysis  capability.  In  addition,  the  specialized  features 
necessary  to  determine  stress  intensity  factors  have  also  been  incorporated  into 
the  code. 

Since  the  actual  16-inch  vessel  was  of  a  jacketed  design,  special  modeling 
methods  were  needed.  Under  elastic  conditions,  the  total  stresses  in  the  liner 
of  the  vessel  were  calculated  as  the  superposition  of  the  residual  compression 
due  to  the  shrink  fit  plus  the  tension  produced  by  the  internal  pressure.  This 
was  modeled  by  solving  two  separate  problems.  First,  the  shrink  fit  was  pro¬ 
duced  by  modeling  the  liner  and  the  external  jackets  as  separate  structures, 
tying  them  together  using  multiple-point  constraints,  and  then  applying  a  ther¬ 
mal  load  that  expanded  the  liner  into  the  jackets.  This  produced  the  shrink  fit 
compression  in  the  liner  in  exactly  the  same  manner  as  the  physical  loading  on 
the  liner.  For  the  internal  pressure  case,  the  vessel  was  modeled  as  above 

1 "ABAQUS  User's  Manual,"  Version  4.5(a),  Hibbitt,  Karlsson,  and  Sorensen,  Inc  , 
Providence,  RI,  July  1985. 
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(liner  and  jacket  separate,  then  tied  together)  and  internal  pressure  was 
applied  at  the  ID  and  along  the  crack  faces. 


RESULTS 

Material  Property  Measurements 

The  fracture  toughness  measurements  are  presented  in  Table  I.  As  is 
readily  apparent,  the  two  liners  have  substantially  different  properties.  The 
reasons  are  explained  by  the  metal  logical  examination  of  these  liners  conducted 
by  Tauscher  (ref  2).  One  liner  had  a  predominantly  ferrite/bainite 
TABLE  I.  FRACTURE  PROPERTIES  OF  LINER  MATERIALS 


Specimen 

Identification 

Test 

Temperature 

("C) 

Fracture 
Toughness 
(Ksi/Tn. ) 

Micostructure 

Valid 

(Y/N) 

A-4 

25 

»  90.4 

F/B 

N 

A-2 

25 

»106.7 

F/B 

N 

B-1 

25 

90.4 

P 

Mar 

B-4 

25 

82.0 

P 

Mar 

B-2 

-40 

63.1 

P 

Y 

B-3 

-40 

58.9 

P 

Y 

A-1 

-40 

108.9 

F/B 

N 

A-3 

-40 

89.7 

F/B 

Mar 

Key; 

P:  Pearlitic  Microstructure 

F/B:  Ferrite/Bainite  Microstructure 
Mar:  Marginal ly  Valid 
1  Ksi/TnT  =  1.099  MPat4 


'^S.  Tauscher,  "Metallurgical  Assessment  of  Two  Liner  Materials,"  Memorandum  For 
the  Record,  Engineering  Support  Division,  Benet  Laboratories,  U.S.  Army  ARDEC, 
Watervliet,  NY,  July  1987. 
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microstructure  (Liner  A),  while  the  second  had  a  predominantly  pearlitic 
microstructure  (Liner  B).  The  vast  difference  in  fracture  properties  was 
expected  since  it  was  well  known  that  the  fracture  properties  of  the  class  of 
steels  examined  with  pearlitic  microstructures  were  very  inferior  to  the  frac¬ 
ture  properties  of  the  same  steel  at  the  same  strength  with  a  ferrite/bainite 
microstructure.  The  fracture  toughness  properties  presented  in  the  table  are 
the  values  obtained  using  the  standard  methods  of  analysis  found  in  Reference  3. 
These  numbers  were  invalid  for  the  more  ductile  material  and  only  marginally 
valid  for  the  more  brittle  case.  The  actual  properties  were  probably  somewhat 
higher  for  the  pearlitic  liner  and  much  higher  for  the  ferrite/bainite  liner. 

The  values  in  Table  I  were  used  to  predict  the  critical  defect  sizes  for  the 
vessels,  since  the  use  of  these  values  yields  a  conservative  prediction. 

Although  the  fracture  properties  were  quite  different,  the  fatigue  crack 
propagation  properties  of  the  two  liners  were  not.  The  results,  presented  in 
graphical  form  in  rigures  2  and  3,  attempted  to  account  for  mean  stress  effects 
or  R  ratio  effects.  This  was  necessary  if  the  effects  of  the  residual  stresses 
from  the  shrink  fit  were  to  be  included  in  predicting  the  remaining  fatigue  life 
of  cracked  liners. 

Figure  2  presents  the  crack  growth  data  in  the  classical  way.  The  crack 
growth  rate  (dl/dN)  is  plotted  against  the  applied  stress  intensity  factor  range 
(AK).  The  straight  line  on  the  log-log  plot  suggests  that  the  material  property 
followed  a  power-law  relationship  and  that  there  was  a  substantial  effect  due  to 
the  R  ratio.  Furthermore,  it  seems  that  there  was  no  effect  on  fatigue  crack 
propagation  due  to  the  different  microstructures.  The  R  ratio  effect  was 

^"Standard  Test  Method  for  Plane-Strain  Fracture  Toughness  of  Metallic 
Materials,  ASTM  E-399,"  Annual  Book  of  ASTM  Standards,  Vol.  03.01,  American 
Society  for  Testing  and  Materials,  Philadelphia,  PA,  1986. 
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expected  because  when  R  is  -1,  substantial  compressive  loading  is  applied  and 
this  tends  to  retard  fatigue  crack  propagation.  In  order  to  account  for  the 
crack  growth  retardation,  the  data  are  plotted  differently  in  Figure  3.  In  this 
case,  the  crack  growth  rate  is  plotted  against  the  maximum  stress  intensity  fac¬ 
tor  applied  during  the  fatigue  cycle  It  is  readily  observed  that  this 

formulation  accounts  for  the  R  ratio  effect  that  was  measured.  All  of  the  data 
fell  in  a  small  scatter  band.  For  crack  propagation  predictions,  an  upper  limit 
power-law  that  represents  the  data  in  Figure  3  is  also  plotted  in  the  figure. 

The  equation  for  this  power-law  is 

dl/dN  =  9.8  X  10-’“  Kn,ax2.86  (1) 

The  unit  of  crack  growth  rate  is  inch/cycle  and  the  unit  of  "is  Ksikin. 

When  the  units  of  crack  growth  rate  are  m/cycle  and  the  units  of  K^ax  3''® 

MPa/in,  Eq.  (1)  becomes 

dl/dN  =  1.9  X  10-”  Kmax2.86  (la) 

Stress  Analysis  Results 

The  K  solutions  developed  from  the  finite  element  analysis  are  plotted  in 
Figure  4  for  the  internal  pressure  case  and  in  Figure  5  for  the  shrink  fit 
loading  case.  In  both  plots,  the  normalized  stress  intensity  factors  are  given 
as  a  function  of  normalized  crack  depth.  The  value  1  is  the  physical  crack 
depth,  and  W  is  the  liner  wall  thickness  (b-a).  The  stress  intensity  factor 
normalization  is  K/a0(iD)ynl,  where  K  is  the  stress  intensity  factor  and  00(10) 
is  the  tangential  stress  at  the  ID  for  the  particular  loading  condition.  For 
the  case  of  internal  pressure  loading,  00(10)  is  the  Lam6  tangential  stress  at 
the  ID  plus  the  pressure  that  acts  on  the  crack  faces  (ref  4) 

^T .  E.  Davidson  and  D.  P,  Kendall,  "The  Design  of  Pressure  Vessels  for  Very  High 
Pressure  Operation,"  The  Mechanical  Behavior  of  Materials  Under  Pressure, 
(H.Ll.D.  Pugh,  ed.),  Elsevier  Publishing  Company,  Amsterdam,  1970. 
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<70(10) 


2pk* 

k2"="I 


(2) 


In  Eq.  (2),  p  is  the  internal  pressure  and  k  is  the  overall  radius  ratio  of  the 
compound  cylinder  (the  liner  plus  all  of  the  jackets),  or  k  =  OD  (outer 
diameter )/ID. 

For  the  shrink  fit  loading  case,  the  value  of  09(10)  is  determined  from  the 
following  equation  (ref  4): 


<70(10) 


-2b‘'  E6  ,(c2-bM(b2-a2), 


(3) 


In  Eq.  (3),  b  is  the  interference  fit  radius,  E  is  Young's  modulus  of  the 
material  (both  liner  and  jacket),  6  is  the  amount  of  radial  interference,  and  c 
is  the  outside  radius  of  the  overall  built-up  cylinder  (OD/2). 

To  facilitate  the  use  of  the  K  solutions  for  the  fatigue  life  predictions, 
it  was  necessary  to  fit  an  algebraic  expression  to  the  numerical  results.  This 
was  easily  accomplished  by  the  use  of  the  normalized  solutions.  We  assume  that 
K  follows  an  equation  as 

K  =  Oq[ID)ViP\  f(l/W)  (4) 

For  internal  pressure  loading, 

f(l/W)  =  1,12  -  0.175  1/W  -  0.515  (l/W)^  +  1.43  (l/W)^  (5) 

For  shrink  fit  loading, 

f(l/W)  =  1.12  -  0.217  1/W  -  0.268  (l/W)^  +  0.943  (1/W)’  (6) 

Equations  (5)  and  (6)  are  plotted  along  with  the  K  solutions  they  represent 
in  Figures  4  and  5.  As  can  be  seen  in  the  figures,  the  agreement  between  the 
equations  and  the  numerical  results  was  excellent. 

^T.  E.  Davidson  and  D.  P,  Kendall,  "The  Design  of  Pressure  Vessels  for  Very  High 
Pressure  Operation,"  The  Mechanical  Behavior  of  Materials  Under  Pressure, 
(H.Ll.D.  Pugh,  ed.),  Elsevier  Publishing  Company,  Amsterdam,  1970. 
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FATIGUE  CRACK  GROWTH  PREDICTIONS 


Using  the  measured  fatigue  and  fracture  properties  along  with  the  K  solu¬ 
tions  generated,  the  remaining  life  of  a  cracked  16-inch  liner  can  be  predicted. 
This  is  done  by  rearranging  Eq.  (1)  to  a  form  in  which  it  can  be  integrated 


AN 


dl/(9.8  X  10-’“ 


*^max2‘86) 


(8) 


In  Eq.  C),  1  .j  is  the  initial  crack  depth.  If  is  the  final  crack  depth,  and  AN 
is  the  number  of  cycles  required  to  propagate  the  crack  from  1^  to  If.  The 
integral  in  Eq.  (7)  cannot  be  integrated  exactly,  but  can  be  easily  evaluated 
numerically  using  Simpson's  rule  or  some  other  technique.  This  was  done  for  the 
specific  loading  of  a  cracked  16-inch  liner. 

There  are  several  specific  variables  that  must  be  fixed  for  the  integration 
to  be  performed.  First,  an  initial  crack  depth  has  to  be  established.  We 
assumed  from  heat  checking  that  a  crack  of  about  0.020  inch  existed  after  the 
vessel  was  cycled  once.  Second,  the  pressure  has  \;o  be  established.  There  were 
two  pressure  ranges  that  could  be  applied  to  these  vessels:  (1)  an  internal 
pressure  of  49.7  Ksi  (342.7  MPa),  and  (2)  an  internal  pressure  of  43  Ksi  (296.5 
MPa).  Third,  the  amount  of  shrink  fit  that  was  in  the  cylinder  has  to  be 
established.  In  the  highest  stressed  region,  the  required  amount  of  inter¬ 
ference  was  dimensioned  as  0.020  +  0.001  inch  (0.51  +  0.03  mm);  0.020  inch  was 
used  as  the  worst  case. 

The  specific  K  solution  for  the  lower  pressure  case  (43  Ksi)  is  given  in 
Figure  6  and  that  for  the  higher  pressure  case  (49.7  Ksi)  in  Figure  7.  These 
figures  demonstrate  how  the  total  K  solution  was  obtained.  The  stress  intensity 
factor  for  internul  pressure  alone  was  determined  for  various  crack  depths.  The 
stress  intensity  factor  for  the  shrink  fit  loading  was  determined  at  the  same 


crack  depths.  The  total  K  solution  was  found  by  point-by-point  addition.  As 
was  expected,  the  49.7  Ksi  internal  pressure  case  was  worse  than  the  43  Ksi 
case.  Using  80  Ksi/Tn.  (88  MPa/iii)  as  the  fracture  toughness  of  liners,  the 
critical  crack  depth  for  a  vessel  subjected  to  49.7  Ksi  pressure  was  0.25  inch 
(6.4  mm),  while  for  the  same  vessel  subjected  to  43  Ksi  pressure,  the  critical 
crack  depth  was  abou':  0.38  inch  (9.7  mm).  If  there  were  no  shrink  fit  (a 
possible  worst  case),  the  critical  crack  depths  for  49.7  Ksi  pressure  would  be 
about  0.13  inch  (3.3  mm)  and  for  43  Ksi  pressure,  about  0.20  inch  (5.1  mm). 

These  were  the  crack  deaths  at  which  the  integration  of  Eq.  (7)  was  terminated. 

The  results  of  the  integration  (the  predicted  crack  growth)  are  shown  in 
Figure  8.  In  this  plot,  three  cases  are  presented:  49.7  Ksi  pressure  and  43  Ksi 
pressure,  both  with  a  0.020-inch  interference,  and  49.7  Ksi  pressure  with  no 
interference.  Again,  these  plots  were  generated  assuming  that  the  initial  crack 
length  was  0.020  inch.  Assuming  that  fracture  occurred  at  the  critical  crack 
deoths  determined  above,  the  remaining  lives  for  the  three  cases  are  as  follows: 
for  49.7  Ksi  with  full  interference  fit,  the  remaining  life  was  3500  cycles;  for 
43  Ksi  with  full  interference  fit,  the  remaining  life  was  6700  cycles;  and  for 
49.7  Ksi  and  no  interference,  the  remaining  life  was  about  1500  cycles. 

The  plots  can  also  be  used  to  determine  the  remaining  life  if  the  initial 
crack  depth  is  more  than  0.020  inch.  For  example,  suppose  that  the  measured 
existing  crack  is  0.050  inch  (1  mm).  We  determine  the  remaining  life  for  the 
49.7  Ksi  case  as  follows:  the  critical  crack  depth  remains  at  0.25  inch  and  the 
number  of  cycles  to  propagate  a  crack  from  0.020  inch  to  0.25  inch  is  3500 
cycles  as  from  above.  The  number  of  cycles  required  to  grow  a  crack  from  0.020 
■inch  tu  0.050  inch  for  this  case  is  about  1500  cycles.  The  1500-cycle  life  has 
already  been  consumed  leaving  the  difference  as  the  remaining  life  for  a  crack 
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to  grow  from  0.050  inch  to  0.25  inch  or  remaining  life  =  3500  -  1500  =  2000 
cycles.  The  same  procedure  can  be  followed  for  the  other  cases  and  any  initial 
crack  depth. 

CONCLUSION 

Fracture  mechanics  has  been  used  to  determine  the  resistance  to  the  pres¬ 
ence  of  cracks  in  1945  vintage- jacketed  thick-walled  pressure  vessels.  The 
procedure  required  the  measurement  of  fracture  toughness  and  fatigue  crack 
growth  rates  for  two  liner  materials  and  the  development  of  stress  intensity 
factors  for  two  loading  conditions.  The  master  crack  growth  rate  curves  devel¬ 
oped  can  be  used  to  assess  the  damage  tolerance  of  these  vessels,  and  to  predict 
the  remaining  life  with  great  accuracy. 


\ 
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Figure  1,  Schematic  of  a  16-inch  pressure  vessel  (1 


a  =  8.0  in. 
b  =  10.16  in. 
c  =  24.5  in. 

k  =  c/a  =  3.06 


in.  =  2.54  cm) . 
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rate  versus  AK  (1  Ksi/in.  =  1.099  MPa/f 
=  2.54  X  lO'i*  m/cycle). 
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Figure  3.  Crack  growth  rate  versus  (1  Ksi/in.  =  1.099  MPa/ni; 

1  Min. /cycle  =  2.54  x  lO'"  m/cycle). 


Figure  4.  Normalized  K  solution  for  internal  pressure  loading. 
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Specific  K  solution  for  43  Ksi  internal 
pressure  (1  Ksi  =  6.89  MPa). 
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